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ABSTRACT 

The disruptive effect of galactic tides is a textbook example of gravitational dynamics. However, 
depending on the shape of the potential, tides can also become fully compressive. When that is 
the case, they might trigger or strengthen the formation of galactic substructures (star clusters, 
tidal dwarf galaxies), instead of destroying them. We perform A-body simulations of interacting 
galaxies to quantify this effect. We demonstrate that tidal compression occurs repeatedly during a 
galaxy merger, independently of the specific choice of parameterization. With a model tailored to 
the Antennae galaxies, we show that the distribution of compressive tides matches the locations and 
timescales of observed substructures. After extending our study to a broad range of parameters, we 
conclude that neither the importance of the compressive tides (ss 15% of the stellar mass) nor their 
duration (~ 10 7 yr) are strongly affected by changes in the progenitors' configurations and orbits. 
Moreover, we show that individual clumps of matter can enter compressive regions several times in 
the course of a simulation. We speculate that this may spawn multiple star formation episodes in 
some star clusters, through e.g., enhanced gas retention. 

Subject headings: galaxies: evolution — galaxies: interactions — galaxies: starburst — galaxies: star 
clusters — stars: formation 



1. INTRODUCTION 

That the violent collision at the origin of a galaxy 
merger acts to develop strong tidal features is a long- 
acknowledged fact of its evolution in time, as en- 
capsulated in the Toomre s equence (see several il- 
lustrations in [Sandage j |196U iToomre fc Toomrej 119721: 
Toomrelll977t IHibbard fc van Gorkomlll996t ILaine et al.l 
20031 lElmegreen et all 120071 ). The burst of star for- 
mation triggered during such an event is a cata- 
lyst for the evolution of entire galactic stellar pop- 
ulations (iBarton. Geller. fc Kenyan! 2000: Bridg eet al.l 
2007i lLTet"al.ll2008fh This" process would be even mor e 
significant at high redshift (z > 1, Madau et al. 1996), 
through either repeated accretion during the formation of 
a major galaxy, or enhanced merger rates driven by the 
formation of cosmological large-scale structures. How- 
ever, the lack of resolution at these redshifts is an intrin- 
sic difficulty when trying to identify stellar associations 
down to a few thousand solar masses, of a size spanning 
a few parsecs. At today's resolution ~ 0.01", a full cov- 
erage of all modes of star formation triggered during a 
merger limits surveys to the low redshift, local Universe. 

Over the past decade, high-resolution HST observa- 
tions of the proto-typical merger NGC 4038/39 (the 
Antennae) have revealed a surprisingly high number 
of dense, massive young clusters of stars in its cen- 
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tion which emerges from such studies pertains to the 
likelihood that stellar associations will survive in the 
rapidly-varying background potential of the merger, once 
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they have outlasted the early phase of internal evolution 
(jElmegreen fc Efremo"vlll997t iBastian fc Goodwinl 120061 : 
iGilbert fc Grahamll2007f). Based on samples of some 10 3 
Antennae clusters. iFall et ail jToO-Yi derived a feature- 
less power- law age distribution (dN/dr oc r _1 ) with 
a median value of ~ 10 7 yr, a trend they argued to 
be independent of mass. Their analysis may be inter- 
preted as implying that young clusters dissolve rapidly 
into their gas cloud nursery, prio r to reaching virial equi- 
librium (IWhitmore et all 120071 : Ide Griis fc Parmentierl 
l2007h . IGielesI (12009) has shown that mass-dependent 
models of cluster evolution leading to dissolution would 
inevitably introduce a feature in the age distribution of 
A ntennae clusters, lending support to the conclusions 
of IFall et alJ (I2005T) . Since cluster dissolution in equi- 
libriu m galaxies invariably lead to Gaussian dis tribu- 
tions (|Vesperinilll997l : iBaumgardt fc Makindl2003h . one 
is drawn to ask whether the power-law properties of 
young cluster distributions in rapidly evolving, out-of- 
equilibrium mergers are shaped mainly by gravitational 
effects, or hydrodynamics, or both. Clearly the evolu- 
tion of the large-scale gravitational potential imprints the 
smaller scales of star formatio n controlled by the hydro- 
dynamics of accretion seeds (|Klessen. Burkert. fc Batd 
119981: iMac Low fc Klessenll2004| ). Knowing how to con- 
nect the two would be a major step toward a global 
understanding of the evolution of stellar populations in 
mergers. 

One way of addressing this problem consists in recre- 
ating observations numerically. Simulations of galaxy 
mergers present much more complicated kinematics that 
increase the complexity level of star formation mod- 
eling when compared, say, to galaxies taken in iso- 
lation. To resolve the many star-forming regions in 
a galaxy merger, several studi es rely on smoothed 
parti cle hydrodynamics (SPH, iGingold fc M onaghan 
Il977h simulations, assuming a star formation recipe 
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mainl y based on a loca l dens it y or pressure t hresh- 
old ((Mihos fc Hernauistl (19961: Bekki fe Chibal 120021 



Bekki et all l2002t iSpringel fc Hernquistj 120031: iLi et al l 
2004t iRobertson fc Kravtsovl 12008 1: iKarl et all 120081: 
Johansson. Naab. fe Burkertj l2009|) . Other a pproaches 



like sticky particles fe .g. iBournaud et al.ll2008f) or adap- 
tive mesh refinement |Kim, Wis e, fc AbeW2009l ) are also 
employed to investigate the formation of substructures 
(star cluste rs or tidal dw arf galaxies, TDGs) at high 
resolution. iBarned (|2004f ) compared the common SPH 
star formation rule with a new prescription based on 
the energy dissipation through shocks, with an appli- 
cation to NGC 4676 (the Mice). He showed that to 
gauge correctly the radiative energy from shocks and to 
recover a more realistic star-formation rate, it is nec- 
essary to take into account the divergence of the back- 
ground velocity field, a strong hint that large-scale mo- 
tion bea rs down on the small-sc ale physics of star for- 
mation. iDobbs fc Pringld (|2009D recently showed that 
the surface density of neither the total gas (H I + H 2 ) 
nor the warm gas only (H I) were good tracers for the 
local star formation rate, considering an isolated disk. 
They also underlined that a strong spiral shock would 
increase the local density, but the velocity dispersion 
too, which preempts a one-to-one relation between star 
formation (averaged over the whole disk) and shocks. 
These results suggest to consider (on top of a den- 
sity contrast) dynamical properties such as the conver- 
gence of the velocity field, either through hydrodynam- 
ical shock dissipation or through a background gravita- 
tional pot ential . On a small er scale, s evera l studies start- 
ing with lEbertl (|1955l ) and iBonnor] (|1956l ). investigated 
the hydrodynamics of gaseous systems like giant molec- 
ular clouds confined by an external pressure (see e.g. 
Kumai fc Tosalll986l: lElmegreenl IrMl Harris fc Pudritd 
1994 lElmegreen fcEfremovl 119971 ) Thev showed that 
such an external effect could trigger the collapse of 
the cloud, so linking galactic-scales hydrodynamics and 
cluster-size phenomena. Despite this impressive body 
of work, the high-resolution needed to link gravitational 
galactic dynamics to the very small scales of star forma- 
tion is still out-of-reach of current 3D numerical models. 

A possible way out of this dilemma would be to im- 
pose the right boundary conditions on a small volume 
within which high-resolution hydro codes may resolve the 
physics of star formation. This paper takes a step in that 
direction and investigates how thes e boundary conditions 
might evolve as a function of time. IRenaud et alj (2008) 
already explored the evolution of the tidal field during 
a merger event, with a model tailored to the Anten- 
nae galaxies. They showed that fully compressive tides 
(see Section [5]) develop during the interaction of two disk 
galaxies and last for few xlO 7 yr. They concluded that 
such a tidal mode could play an important role in trig- 
gering the formation of star clusters, or at least prevent 
them from destruction. Here, we detail further and ex- 
tend these results through a larger sample of progenitor 
galaxies and merger orbits. It is hoped that this more 
general yet precise description of the tidal field will in 
turn allow more realistic boundary conditions for vol- 
umes where star formation takes place. In Section 2, we 
present analytically how compressive tides may develop 
when two potentials overlap, as in the case of a major 



galaxy merger. Section 3 describes the application to 
A^-body simulations and the limitations of the method. 
Section 4 details the Antennae reference model and Sec- 
tion 5 makes direct comparisons between the parameter 
survey and this reference. We close with a discussion on 
the link between the physical process described and the 
observational facts. 

2. GRAVITATIONAL TIDAL TENSOR 

The tidal field describes how the internal dynamics of 
an extended body, e.g. a star cluster, is stirred by the 
influence of an external source of gravitation, say the 
host galaxy. In this case, the differential force due to the 
galaxy which exists between any two points within the 
cluster, and separated by £ is given roughly by SF = IT 
where T is the tidal tensor of components 



T» EE -di (dj 



(1) 



which are minus the second derivatives of the gravita- 
tional potential <p along the j-th and «-th spatial axes. 
Because the tidal tensor is symmetric and real-valued, 
it can be set in orthogonal form. In this representation, 
the eigenvalues A gauge the strength of the tide along the 
corresponding eigenvectors. According to the sign of an 
eigenvalue, one may characterize the tide down the asso- 
ciated eigenvector to be compressive (A < 0) or extensive 
(A>0). 
The trace of T reads 



J2 -d> {di<t>) = -V 2 = -iirGp, 



(2) 



where p is the density and the last steps follows from 
Poisson's law. Since the density is a positive-definite 
quantity at every point in space, it is impossible to com- 
pute three positive eigenvalues simultaneously. Hence, 
the number of axes along which the eigenvalue may be 
positive ranges from zero to two. Cases for which one 
or two eigenvalues are positive occur in well-known con- 
figurations (Earth-Moon inte raction, galac t ic or cluster 
tides , spiral waves, see e.g. IValluril [1991 : IDekel et al.l 
120031: IGieles et all 120071) . However, the case where all 
three eigenvalues are negative has been somewhat over- 
looked in the literature. When that occurs, the tidal 
force field converges to a common point from all direc- 
tions, and so the tides are fully compressive with respect 
to that point. Hereafter, we will refer to such a situation 
as a compressive mode of the tidal field. 

From a mathematical standpoint, one may remark that 
a compressive mode only exists in the convex intervals of 
the twice-differentiable gravitational potential 4 . A sim- 
ple analytical e xercise demonst rates indeed that co red 
potentials, e.g. IPlummeri (|1911l) and IDehnenl (|1993L for 
< 7 < 1), show a fully compressive region near their 
cente r of mass, while p rofiles with a cusp, e.g. point - 
mass. lHernquisti (fl9 90). Dchncn (1993, for 1 < 7 < 3) or 
NFW (|Navarro et al.|[l997f ). have extensive tides along 
at least one axis everywhere in space. (Table Q] lists the 
analytical expressions of the cases we considered.) 

The overlay of two identical density profiles (cuspy or 
not), one shifted with respect to the other, can modify 

4 The tidal tensor is the Hessian matrix of the potential and 
thus, vanishes at the inflection points. 
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TABLE 1 

Tidal tensors of typical density profiles 



Profile 



(x = r/ro) 



T xx {x = r/r ,y = 0,z = 0) 



Point mass 



Plummcr 



Hcrnquist 



Dehnen 



NFW 



GM 
x 

GM 



to yl + x' 2 

GM 
ro(l + x) 



-GM 



2GM 

x 3 

1 - 2x 2 

(1 + x 2 ) 5/2 

2GM 



GM 



ro(2 - 7) 
Gp 







1 - 


(l + a:) 







GM 



4(l + xf 

2x + 7 - 1 
r$ x~* (1 + x) A —l 



■ In (1 + x) 



Gp 



21n(l + x)- 



3x 2 + 2a; 



(i + xY 



* For 7 = 2, the potential is given by GM/ro ln[x/(l + x)] but the expression 
of the tidal tensor remains the same. 



the potential well such that compressive modes now ap- 
pear where none was present in isolation. In Fig. [TJ we 
show the examples of Plummer (top) and Hernquist (bot- 
tom) potentials, in isolation (left) and when two identi- 
cal profiles overlap, one shifted along the vertical axis 
(right). The color code on that figures gives blue for ar- 
eas where at least one tidal component is extensive, and 
red for a fully compressive mode. The Plummer potential 
in isolation has a large volume about the center where 
the tidal field shows a compressive mode; this situation 
does not occur in the case of a Hernquist profile. By 
contrast, compressive modes exist in both cases where 
two potentials overlap (right-hand frames). This demon- 
strates that even for the case when an isolated galaxy 
hosts no compressive tides anywhere, such tides may still 
develop widely as two galaxies come closer and merge. 

The situation depicted in Fig. [T] will occur naturally 
in the course of a galaxy merger. When a progenitor 
galaxy (sum of a disk, bulge and dark matter halo) comes 
close enough to a second, similar galaxy, local cores in 
the global potential may appear and lead to compressive 
modes of tide that may last for a significant period of 
time. As a side issue, it follows that a compound model 
of a progenitor galaxy may itself display regions where 
tides are compressive although each component, taken 
individually, may not. This phenomenon would in turn 
affect the matter (stars, clusters, interstellar medium) 
inside these areas and possibly alter their dynamical and 
chemical properties (e.g., through enhanced retention of 
metals). However, the very existence of such tidal modes 
will be impossible to predict in a realistic simulation of 
a galaxy merger with analytic tools only, owing to the 
irregular morphology of the system. 

3. NUMERICAL METHOD 

While simple density profiles allow an analytical ex- 
amination of the tidal field, the same approach applied 
to mergers is clearly out of reach. To overcome this, we 
set up three-dimensional ./V-body simulations to explore 
the time-evolution of the gravitational potential and its 
associated tidal field. This section gives details of the 
numerics used in this work. 

3.1. Method 






compressiv 


6 Jr. 


extensive 



-1.0 -0.5 0.0 0.5 1.0 

Fig. 1. — Maximum value of the tidal tensor. The Plummer po- 
tential (top) shows an intrinsic compressive region (red) in isolation 
(top left) as well as for an overlap of two similar potentials (top 
right) separated by 1.5 ro (where ro is the core radius). For the 
Hernquist potential (bottom), the overlap of the two wells (bottom 
right) separated by 0.4 ro also presents a compressive region. How- 
ever the isolated case (bottom left) is completely extensive. The 
dashed circles represent the characteristic radii ro- 

For simplicity, we used purely gravitational models (no 
gas). This approach holds good as long as the gas rep- 
resents a small fraction of the total mass of the galaxy, 
which is usually the case in the local Universe. Therefore 
neglecting the gas component should not affect signifi- 
cantly the shape of the gravitational potential and the 
characteristics of its tidal tensor. The numerical integra- 
tion of orbits was performed with the public version of 
Walter Dehnen's code gyrf alcDN with a compact kernel 
K0 (see lDehnerJl200l) . The mass- and length resolutions 
were linked to typical globular cluster scales (i.e. 10 5 M 
and 40 pc). All the particles have the same mass to re- 
duce two-body relaxation. 
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In our compound A^-body galaxy models (disk + bulge 
+ dark matter halo), we focused on disk particles when 
computing the tidal tensor, since, being dynamically 
cold, they give rise to large-scale arms, bridges, and other 
transient features. To compute the tensor T, first a mass- 
less particle is placed on each of the six faces of a cube 
centered on a given disk-particle coordinates. In a sec- 
ond step, the force field is computed at the po sition of 
each pseudo-particle using the getgravity tool (|Dehnenl 
2002). Note that the gravity of the disk particle sit- 
ting at the centre of the cube is subtracted explicitly. 
Finally, we implemented a first-order-accurate differenc- 
ing scheme to compute the gradients of the forces by 
pairing the pseudo-particles as needed to derive the six 
independent components of the symmetric tidal tensor. 
The eigenvalues and eigenvect ors we r e then computed 
efficiently using the method of iKoppI (|2008l ). The sign 
of the largest eigenvalue determines whether the central 
disk particle sits in a fully compressive volume, or not. 

In order to obtain statistically significant results, we 
applied the above procedure to all the disk particles of 
every simulation snapshot (each separated by a time in- 
terval of ~ 2.5 Myr). In that fashion we retrieved the 
evolution of the maximum eigenvalue to trace the his- 
tory of the tidal field of disk particles, a binary test on 
the character of the tide (compressive or extensive) expe- 
rienced by a particle throughout the merger simulation. 

3.2. Tests 

Using the analytical Plummer model (see Table [1]) as 
a test case, we computed the gradients of the force (i.e. 
the tidal tensor) using a first-order-accurate finite dif- 
ference scheme. The size of the cube 8 defines a length 
for the differentiation scheme which should be chosen to 
minimize errors. Large values of 8 will fail to resolve 
fine features in the potential, while very small values are 
prohibited when round-off errors kick in. In practice, 
however, errors attributable to the difference scheme ac- 
counted for less than 15% of the total uncertainties: they 
have been omitted from the remainder of the discussion, 
for the sake of clarity. 

The main source of errors comes from the TV-body ren- 
dition of the potential. When integrating the orbits of an 
A^-body model, a softening length e is often invoked to 
suppress small-scale errors in the force field. This length 
also sets a resolution scale for all quantities bearing on 
the dynamics. Thus, one must choose S > e to compute 
the tidal tensor coherently with the resolution of the sim- 
ulation. The more tricky issue of finding a suitable upper 
bound for d for a given e remains. 

To address this question, a series of tests were per- 
formed with an TV = 5 x 10 5 bodies realization of a Plum- 
mer sphere. Eigenvalues and eigenvectors were com- 
puted throughout the model with the numerical scheme 
outlined above and for various values of 8. The best 
results (Fig. ^ compared to the analytical values had 
8 m 0.08 ro « 0.05, where ro is the Plummer core radius. 

This value of 8 may be compared to the softening e 
of the force field (|Merrittill99l lDehnenll2001h which has 
been set to 0.01 for the choice of N adopted here. There- 
fore, an empirical relation suggests itself 8 « 5e. 

To determine the scope of this relation, we repeated 
the calculations outlined here to several compound mod- 
els of the Antennae galaxies, marking a clear break with 
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Fig. 2. — Maximum eigenvalue of the tidal tensor of a Plummer 
sphere along one radius. A cube size of 8/ro = 0.08 (with ro, the 
core radius) minimizes the dispersion (<r ~ 0.13) relative to the 
analytical solution (see Table [TJ , while smaller and larger values 
lead to less precise fits (e.g. a size S/ro = 0.02 shows a ~ 0.19). 
Small cubes increase the noise while larger ones do not match the 
analytical solution. For a Plummer distribution, the compressive 
region is a sphere of radius rrj/v2, which encloses ~ 19% of the 
total mass. 

the spherically symmetric Plummer potential, with val- 
ues of N ranging from 350,000 to 1.4 x 10 6 . We found 
that the relation 8 — 5e holds well over this range of par- 
ticle number. It was therefore applied to all simulations 
presented in the following parameter survey. 

4. THE ANTENNAE AS A REFERENCE 

The starting point of our su rvey is the model of the 
Antennae galaxies presented in Rcnaud et al.l (|2008l ). It 



(Boilv et al. 


2001) 


Hernauisti ( 


1990) 



which gathers an exponential disk, a 
bulge and an isothermal dark matter halo. The mass 
ratios and length-scales of each component have been set 
to get two SO galaxies. Initial parameters are described 
in Tabled 

A total of 1,400,000 particles have been used. As men- 
tioned before, our resolution matches the typical cluster 
scale (e ~ 40 pc and 10 5 M Q ). The progenitor galax- 
ies have been set on a bound orbit with an eccentricity 
e w 0.96 and an initial separation of 75 kpc. Accord- 
ing to the relation 8 = 5e (see Section 13.21) , we set the 
size of our cubes to 220 pc, i.e. well beyond the typical 
half-mass radius of a star cluster. 



4.1. Time evolution 

Fig. [3] tells the history of the interaction. Most of the 
tidal features originate from the disks during the first 
pericenter passage. The two long tails expand as the 
galaxies move apart from each other. The central regions 
are strongly mixed during the second passage and brake 
significantly, which leads to an oscillatory phase and fi- 
nally a rapid merger. The best match with observations 
(morphology and kinematics) occurs at the beginning of 
the second passage (Fig. [31 panel h), some 300 Myr after 
the first one. 

The evolution of the compressive tidal mode (in term 
of intensity and mass fraction) clearly shows the same 
steps. Fully compressive regions are displayed as red 
dots in Fig. [3j Initially, they are located in the intrinsi- 
cally compressive regions at the very center of the galax- 
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TABLE 2 

Parameters of the reference Antennae model 



Parameter 


NGC 4038 


in v ' v_. 4uoy 


Numbers of particles 




iV disk 


2 x 10 


2 x 10 5 


AT 

1 * bulge 


1 X 1U 


1 x 10 5 


AT 

iV halo 


4 X 1U 


4 x 10° 


— T' 1 1 i|. 

bcale- lengths 






h 

'•-disk, vertical 




0.2 


^disk 


1.0 


1.0 


f'bulgc 


1.0 


1.0 


Hialo 


7.0 


7.0 


(Jut-oii radii 






^disk. vertical 


2.0 


2.0 


'--disk 


5.0 


3.0 


^bulgc 


1.0 


1.0 


•-halo 


7.0 


7.0 


Mass ratios 






W-disk 


1.0 


1.0 


^bulgc 


1.0 


1.0 


m halo 


5.0 


5.0 


Toomre paramet 


ers 




Q 


1.5 


1.5 


Disks inclinations 




Ox 


60° 


-60° 


By 


0° 


30° 


initial coordinates 






(6.0,6.0,0.0) 


(-6.0,-6.0,0.0) 


{v x ,Vy,V Z ) (" 


0.5, -0.25,0.0) 


(0.5,0.25,0.0) 



Note. — The parameters are given in numer- 
ical units (G = 1) with the conversion factors 
4.4 kpc, 1.9 X 10 2 km.s' 1 and 3.6 X 10 10 M for 
the length, velocity and mass respectively. The non- 
dimensionalization is based on the spatial extension 
of the tidal t ails and the peak in the velocity field 
measured by Hibbard ct al. (2001), adopting a dis- 
tance to the Antennae of 19.2 Mpc. The Keplcrian 
equivalent orbit has an eccentricity of 0.96, a distance 
to pericenter of 1.31 (= 5.78 kpc) and a semi major 
axis of 17.74 (= 78.06 kpc). 
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Fig. 3. — Morphology of the disks of the Antennae model every 
50 Myr. On each row, the black line represents 10 kpc. t = 
(b) corresponds to the first pericenter. The tidal bridges and tails 
are created in the first 10 7 yr after the first passage, while the 
galaxies move apart. The best match with observations is obtained 
at t Si 300 Myr (h). Red dots denote the particles standing in fully 
compressive tidal mode. 
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Fig. 4. — Evolution of the mass fraction of the disks of the Anten- 
nae in compressive mode. The peaks denote the main steps of the 
merger such as the pericenter passages (at t = and t = 300 Myr). 
Lower resolutions give similar results (with deviations from the 
1,400,000 resolution of ct = 0.79 and cr = 1.0 for 700,000 and 
350,000 particles respectively). 

ies (panel a) but, as the potentials overlap, they quickly 
propagate to the tidal arms and bridges (panels c, d, e). 
At t w 150 Myr (panel e), a compressive spiral structure 
is visible in the central part of the disks. Note on this 
panel a large (~ 20 kpc) compressive zone exists near 
the tip of the southern tail, matching the positions of 
the tidal dwarf galaxy candidates identified by several 
authors (|Schweizerl[l97l iMirabel et alJll992ft. We also 
note t hat age estimates of ~ 120 Myr from lHibbard et ahl 
(2005) correspond well to our results. 

During the second passage and the merger phase (pan- 
els g, h, i), the fast evolution of the mass distribution 
mixes the material in the central parts, which creates 
patterns of compressive mode, especially in the overlap 
region and around the two nuclei. Using t his configu- 
ration of the potential, iRenaud et al.l (|2008l Fig. 2) de- 
rived a map of the compressive regions. They used a 
3D grid-based method applied on the potential given by 
the iV-body simulation to get rid of projection effects. 
The good match be tween this map and the candidate 
clusters identified bv lMengel et a l. (2005), in addition to 
the TDGs as discussed above, is a strong hint that the 
compressive tides are to be linked to the formation of 
substructures. 

Fig. H] traces the evolution of the mass fraction in com- 
pressive mode along time. Results for lower resolutions 
(factors 1/2 and 1/4) are also plotted to demonstrate 
that the derived quantities are independent of the num- 
ber of particles used to represent the galaxies. The main 
peaks correspond to the two pericenter passages and 
the smaller spikes are linked to the formation of smaller 
structures such as bridges. 

In the background of all these steps, an intrinsic subset 
(~ 3%) of the stellar mass stays in compressive mode in 
the galaxies' centers (within a radius of ~ 1.3 kpc, when 
the disks are at rest). This offset is due to the initial 
mass distribution and thus, will change with e.g. the 
dark matter halo shape. Because of orbital mixing, some 
particles will stay in this region for a very long time, 
but others will be replaced, keeping this fraction almost 
exactly constant (deviation of 0.2% over 12 Gyr). 

4.2. Characteristic times 
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For an isolated progenitor, the most bound particles 
stand in the central part of the potential well and thus, 
in the intrinsic compressive region. They will stay in this 
region for an arbitrary long time provided their energy 
is smaller than the energy level at the boundary of the 
compressive region. The fraction in this regime (con- 
stantly compressive) represents ~ 1.1% of the total disk 
mass 5 . In the following, we subtracted these objects from 
the statistics, and concentrated on the particles entering 
compressive modes. The others ~ 2% in the compressive 
region (which gathers 2.0% + 1.1% ~ 3.0%) enter and 
leave the central part because of the larger angular mo- 
mentum of their highly eccentric orbits. Therefore, the 
duration of the compressive mode for individual parti- 
cles (in an isolated galaxy) can be studied by considering 
their energy and angular momentum and deriving how 
much time they will spend inside the central compres- 
sive region. A good knowledge of these statistics allows 
to characterize the effect of the intrinsic tidal field, re- 
garding to the influence of the merger. 

At this point, we introduce two timescales. The longest 
uninterrupted sequence (lus) represents the duration of 
the longest continuously compressive episode, while the 
total time (tt) is the count of all compressive snapshots 
which might be separated by periods of extensive tides. 
As a particle would probably undergo a series of on/off 
compressive episodes, the tt should be larger than, or 
equal to, the lus. The limit is reached for the 1.1% con- 
stantly compressive particles whose lus equals tt. A sec- 
ond family of particles counts those which orbits allow 
to enter the compressive region for a short timescale. 
These particles present a medium energy, a low angular 
momentum and thus an intermediate valued lus. The 
last part gathers the high energy and high angular mo- 
mentum particles that merely never enter into the com- 
pressive region. Their lus should be extremely small. 

The value of the angular momentum at the edge of 
the compressive region sets the upper limit of the par- 
ticles which could experience compressive tides, during 
at least one snapshot (2.5 Myr). Particles with a higher 
angular momentum would not enter the central region, 
while those with a smaller value could get into it. The 
last represents 7.7% of the disk mass. Therefore, because 
of orbital mixing, some of these ~ 7.7% of particles will 
temporary join the constant fraction of ~ 1.1% to build 
the measured ~ 3% compressive fraction. To estimate 
the duration of their compressive episodes, we compute 
the velocity associated to the angular momentum at the 
boundary of the compressive zone and derive a mean 
orbital period of 53 Myr. This means that the parti- 
cles enter and leave the central compressive region on a 
timescale of ~ 50 Myr. (Note that the free-fall time of 
this region is iff = 26.3 Myr and that 2tg, i.e. the time 
needed to get to the center and exit the region, is approx- 
imately the value we derived.) Therefore, one should not 
expect a high fraction (of the non-constantly compressive 
particles) to present a lus greater than 50 Myr. 

Fig. [5] reflects the points described above (for lus only). 
It plots the integrated mass fraction M(> At) which 
experienced a continuous compressive event longer than 
At, for the Antennae case as well as for an isolated pro- 

5 Some ~ 1% of these particles will leave the compressive region, 
over 1 Gyr, because of two-body relaxation. 



1 00 00 


\ 


I - 1,400,000 particles 






700,000 particles - 

L 350,000 particles ■ 


10.00 




Isolated 




1 .00 






0.10 












0.01 







50 100 150 200 

Compressive mode interval At [Myr] 

Fig. 5. — Distributions of the longest uninterrupted sequence 
in compressive mode. Integration has been performed for 
t £ [—100,1000] Myr for the Antennae and for the same time 
(1100 Myr) for the isolated case. Again, lower resolutions presents 
the same results (a = 0.48 and a = 0.43 for 700,000 and 350,000 
particles respectively). 

genitor. These distributions are normalized by consider- 
ing the non-constantly compressive particles only, which 
means that M(At = 0) = 100%. As expected, in the iso- 
lated case, ~ 8% of the mass experience a compressive 
episode longer than 2.5 Myr while only 0.1% will be in 
this mode for more than 50 Myr. 

For the Antennae case, the distribution is well fitted by 
exponential laws M(> At) oc exp (At/r) in the time in- 
tervals ATi =]0, 50] Myr and AT 2 = [50, 100] Myr. For 
the lus, we obtain t\ ~ 8.0 Myr and r 2 ~ 50.2 Myr 
for the AT\ and AT 2 , respectively. The isolated case 
presents similar slopes. 

One has to be careful before jumping to conclusions. 
These distributions depend on the time of integration, 
i.e. the time of calculation, during which one allows the 
particles to enter into the compressive regime. For an 
integration over a long time, the distributions would be 
flatter, because a higher number of particles could expe- 
rience one or more compressive episodes. To check the 
validity of our results, we plot these distributions for var- 
ious integration times (from 600 Myr up to 2.1 Gyr) in 
Fig. O for the Antennae case. In addition to the lus, the 
total time (tt) spent in compressive mode (i.e. the sum of 
all the compressive events) is also shown. The position of 
the knee (at ~ 50 Myr) does not vary but, as expected, 
we note a large difference in the tt distribution (up to a 
factor 6). As the purpose of this work is to concentrate on 
the first 100 Myr of these modes, this variation does not 
strongly affect the conclusions presented here. For most 
of the particles, the lus begins at the first pericenter pas- 
sage and thus, the associated distribution remains stable 
to a change of integration time, provided it includes the 
two passages (peaks in Fig. HJ). It is aiS0 important to 
note that a tt distribution integrated over a long period is 
physically less relevant as it is punctuated by longer and 
longer extensive episodes, which could blur the effects of 
the previous compressive events. 

For the Antennae, one could explain the short peri- 
ods with the existence of transient features and link the 
longer episodes to the intrinsic compressive areas, like 
in the isolated case. Fig. [7] confirms this assumption 
by showing the distribution of the compressive periods 
for various regions of the merger. It distinguishes three 
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Fig. 6. — Distribution of the longest uninterrupted sequence (lus) 
and total time (it) spent in compressive regime, for various inte- 
gration times. If a strong dependence is visible for long episodes, 
the short events (< 50 Myr) remain well-defined, whatever the in- 
tegration time is. We note that the position of the knee is not 
affected. 
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Fig. 7. — Distribution of time intervals spent in compressive mode 
(lus only) for different regions of the Antennae galaxies. The tran- 
sient features (tails and bridges) present very short periods (always 
less than 50 Myr), while the nuclei have a long period regime, too. 
Thanks to the symmetry of our Antennae model, the structures of 
the southern galaxy only are presented. 

main areas: the nuclei, the tidal tails and the bridges. 
As expected, the transient structures like the tails and 
the bridges show very short compressive periods while 
the central parts (nuclei) present the knee at ~ 50 Myr, 
between short and longer periods. Therefore, the intrin- 
sic compressive regions are the only ones to experience 
compressive episodes of few 100 Myr. A shorter period 
is also visible and matches the associated distributions 
for the isolated progenitors. Hence, one can tell apart 
the effect the orbital mixing (short and long events) and 
the role of the merger (short periods only). 

4.3. Compressive history 

The typical period of 50 Myr reported above matches 
the time scale of star cluster formation. In addition, the 
positions of the compressive regions matches amazingly 
the location of observed young structures. These clues 
underline the role of the tides in the birth of clusters or 
TDGs. However, one may ask how the statistics pre- 
sented above reflect the actual tidal history of a single 
object. Indeed, the lus concept (introduced for statisti- 
cal purposes) only tells us about one compressive event, 



while a star cluster orbiting in the Antennae is expected 
to experience many of them. Therefore, we propose a 
new description of the local tidal field, by selecting a few 
particles in each region of interest of our Antennae model: 
the central part, the bridges, the tail in the vicinity of the 
nuclei and the TDGs region. The orbits of these particles 
are plotted in Fig. [5] 

This new approach allows to retrieve the tidal history 
of individual orbits and, with the statistics obtained be- 
fore, to get a broader picture of how compressive modes 
evolve. In other words, Fig. [5] shows some realizations 
of the statistics we presented. These examples define 
realistic boundary conditions that will be used for the 
simulation of star clusters. Note that each particle rep- 
resents a special case among the mean behavior of the 
associated region. Therefore, one can pile up snapshots 
during which a body stands in compressive mode, and 
consider that as a compressive event, even if it is punc- 
tuated by short extensive episodes. That is, most of the 
compressive occurrences follow the global history of the 
merger, as related on Fig. [3] and [U 

Indeed, the particles spotted in the bridges (B, C) or 
the tails (F, G) show their major compressive event after 
the first pericenter. One can also track the fall-back of 
the particles into the nucleus (E) or simpler, the long 
period of a particle that stays in the center for the entire 
simulation (A). 

This figure helps to distinguish the major families of 
compressive histories one could find in a merger. The 
behavior of these selected particles matches the statistics 
derived above, in terms of position in space and time, and 
time scale. 

5. PARAMETER SURVEY 

We performed a series of simulations, each with a phys- 
ical parameter set to a different value compared to the 
reference Antennae run, with a view to quantify the im- 
pact of these parameters on the statistics of the tidal 
field. This section describes the various models; Table [3] 
sums up the main results, and Fig. [5] shows the mor- 
phologies and the spatial distribution of the compressive 
regions after the first pericenter passage. 

5.1. Prograde - retrograde encounter 

It has been clear since iToomre fc Toomrd ()1972l ) that 
the large-scale morphology is linked to the initial param- 
eters of the stellar disks. Indeed they showed that, e.g., 
the intrinsic angular momentum of a galaxy (spin) could 
be the key factor in the creation of the tidal tails. We 
propose a similar exploration of the influence of the spin, 
for disks embedded in a dark matter halo. 

The reference Antennae calculation is a prograde- 
prograde (PP) encounter, i.e. both disks' spins are cou- 
pled with the orbital angular momentum of the galaxies. 
This configuration allows the growth of the two observed 
tails. In this model, ~ 30% of a single disk mass is ejected 
to form the associated tidal tai l. (This mas s loss was 
noted in earlier studies, see e.g. IBarnesl lT988.) Varying 
the spin of one or both of the disks allows to change the 
tidal structures created during the encounter. For a ret- 
rograde galaxy (spin anti-aligned with the orbital angular 
momentum) , the bound mass remains approximately the 
same. This affects the ensuing morphology, but also the 
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Fig. 8. — Orbits of the characteristic particles (in the orbital 
plane) for the inner regions (top), intermediate zone (center) and 
large scale (bottom). Dots show the compressive events. The two 
pericenters are marked with crosses (+ and X respectively) on every 
orbit. Particles A and B stay in the vicinity of the central region. 
C goes into one of the bridges, D falls back in the merger remnants 
in the merger phase (t = 300 Myr). Particle E is ejected to the 
southern tail by the first encounter but falls back into the nucleus 
before the second. F stays in the tail and G is part of the TDGs 
region. One position (i.e. one snapshot) corresponds to 2.5 Myr. 
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a Mass fraction in compressive mode added at the first pericenter 
passage (compared to the isolation stage). 

b Mass fraction in compressive mode added at the merger phase. 

c Exponential ti mesc ale of the lus distribution for the time interval 

ATi (see Section IO). 



dynamical mass. Indeed, in a retrograde-retrograde con- 
figuration (RR), the motion of the center of mass of each 
galaxy shows a larger separation (~ 4 kpc) than for PP, 
and thus an orbital period ~ 100 Myr longer. 

We first set a prograde-retrograde (PR) simulation by 
flipping one of the two disks. Only the prograde galaxy 
creates an arm while the mass of the other disk remains 
almost constant during the encounter (see the second 
row of Fig. [9]). At the second passage (i.e. just before 
merger), the amount of material in that disk which can 
experience a compressive mode is much higher. 

A merger event is violent and affects entirely what is 
left of the disk, after the first passage. In a first approx- 
imation, one can consider that the reaction of the disk 
material would be the same, whatever orientation the 
spin takes. However, for prograde galaxies, a significant 
fraction of this mass has been ejected in the form of tidal 
tails. If a given fraction X of the bound mass undergoes 
a compressive mode, the equivalent fraction for the to- 
tal mass (disk + tails) must be lower. In other words, 
a retrograde disk provides more mass to enter in a com- 
pressive mode, thus the equivalent fraction is higher than 
for a prograde configuration. 

For the Antennae run (PP), 14% of the stellar mass 
is in compressive mode at the second pericenter passage 
(t = 300 Myr). Each disk has lost 30% of its initial mass 
into tails. Therefore, X = 0.14/(1.0 - 0.3) = 20% of the 
bound mass stands in compressive mode. For a retro- 
grade galaxy, we expect ~ 20% of the mass to experience 
compressive tides. In the PR case, we have one prograde 
disk (14%) and one retrograde progenitor (20%) of equal 
mass. Therefore we expect (20% + 14%)/2 = 17% of 
the mass in compressive mode. Because of the change 
in dynamical mass, the second pericenter passage of PR 
model should occur ~ 50 Myr after those of the PP con- 
figuration. 

Using the same arguments, we can predict that the RR 
case would have up to 20% of its mass in compressive 
mode (because no tail has been created), and this, ~ 
100 Myr after the PP case. 
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Fig. 9. — Morphology of the 12 models at (approximately) the same dynamical stage, between the first two passages. The same line of 
sight has been used to ease the comparisons. The images are centered on the center of mass of the entire system (dark matter included). 
On each image, the black line corresponds to 10 kpc. Red dots denote the particles standing in fully compressive mode. 



Fig. [10] shows the results. The initial ~ 3% offset 
and the first peak correspond to non-altered disks and 
thus, arc common to the three configurations. When 
the two galaxies move apart (t £ [0, 250] Myr) the de- 
creasing rate of the mass in compressive mode depends 
on the geometry and remaining mass in the central re- 
gions. The second passage, as expected, spreads from 
t = 300 Myr for the Antennae case to t = 400 Myr for 
the RR model. The maximum compressive fraction is in 
good agreement with the argument discussed above. The 
last small peaks stem from the oscillations of the distance 



between the two nuclei which finally merge into a single 
ellipsoid. The three configurations show these spikes, as 
the final merger phase does not depend on the external 
tails. Once the overall potential well stabilizes, the com- 
pressive fraction remains constant, within fluctuations of 
similar amplitude than the initial noise. 

Even if the compressive peaks are shifted in time and 
have different amplitudes, the overall behavior of these 
curves remains the same, whatever the orientation of the 
spins is. However, it is important to remark that in all 
these cases, the impact angle between the two progenitors 
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Fig. 10. — Evolution of the mass fraction in compressive mode. 
Relative amplitude of the peaks and time shifting can easily be 
explained by mass loss consideration. See text for more details. 
Note that the retrograde-prograde configuration (not shown here) 
almost matches the prograde-retrograde one. as our progenitors are 
similar. 
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Fig. 11. — Distribution of the lus in compressive mode for the PP 
(Antennae), PR and RR cases. For the thin curves, the integration 
time has been set from 100 Myr before the first passage, to three 
times the pseudo-period (i.e. the time-lapse between the main two 
peaks). In all cases, the short periods and the knee at 50 Myr are 
similar. For the thick curves, integration has been stopped just 
after the merger event (at t = 375, 440 and 500 Myr for the PP, 
PR and RR cases respectively). 

is the same (modulo 180°). We also run simulations with 
other angles to check this. Even a rotation of 90° of one 
of the progenitors does not affect the mass fraction in 
compressive mode. Only the time interval leading up to 
the second peak increases by 50 Myr. 

A visual inspection of the width of the peaks in Fig. [TO] 
allows to predict that the short time scale in compressive 
mode (lus distribution) would not change with the spin of 
the galaxies. Indeed, Fig.fTTIshows that the slope derived 
by exponential fits in section 14.21 remains the same, as 
well as the position of the knee at ~ 50 Myr. 

Analyzing the longer periods is more difficult, as they 
strongly depend on the integration time. Because of the 
mass loss of the prograde galaxies, the dynamics of the 
three systems are not the same. We define a pseudo- 
period as the time-lapse between two major peaks of 
compressive mass fraction (see Fig. 1101) . For the thin 
curves of Fig. [TTl we set arbitrarily the integration period 
to three pseuso-periods, i.e. 1000, 1150 and 1300 Myr for 
the PP, PR and RR cases. The difference in the final 
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High eccentricity is required to avoid a cen- 
tral collision (due to orbital decay) that would 
hide the effect of the initially finite pericenter 
distance. 

mass fractions visible in Fig. [TO] is integrated for a long 
time. The RR configuration, which shows only ~ 2% of 
the mass in compressive mode after the merger phase, 
has a lower lus distribution than the PP one (which 
equivalent value is ~ 3%). 

However, in order to get a fair comparison between 
the models, one can link the integration time to a 
dynamical stage, and not a period based on the or- 
bit. By stopping the integration just after the merger 
(Fig. [IT] thick curves), we deal with three objects show- 
ing the same dynamical stage in the Toomre sequence 
(|Toomre fc Toomrd [l97l . In that case, the compar- 
isons regarding time are much more relevant. Indeed, 
the lus distributions are very similar, even for long pe- 
riods, showing that the spin of the progenitors does not 
influence the timescales of the compressive tides (see Ta- 
ble©. 

5.2. Orbits 

The next part of this survey focuses on the influence 
of the progenitor orbits, keeping all the other parameters 
(e.g. the spin or the mass ratio) as in the Antennae 
reference. Because the galaxies are not point-masses, it is 
not possible to use a Keplerian description of their orbit. 
However, before the first pericenter passage, the motion 
of the centers of mass of the disks can be approximated 
as a two-body problem. 

To ease the comparisons, in the following we will refer 
to the equivalent Keplerian orbits. Tabled] describes the 
parameters (eccentricity e and pericenter distance d) of a 
set of representative orbits. Because of the orbital decay 
due to dynamical friction, this two-body approximation 
is no longer valid after the first passage. 

The evolution of the mass fraction in compressive mode 
is given in Fig. 1121 with the Antennae for comparison. 

At the first passage, the OP run shows a similar behav- 
ior to the reference. However, its large eccentricity gives 
it a long orbital period, which leads to a second pas- 
sage 6 delayed by ~ 600 Myr. A third encounter occurs 
just before the merger phase and its oscillatory mode. A 
contrario, a small orbit like OE brings back the galax- 
ies very quickly and thus induces the merger while the 
transient tidal features still exist. 

For a distant passage (OD), the overlap of the poten- 
tials which leads to a compressive zone is much smaller 
and thus, the corresponding mass fraction is lower (~ 6% 
instead of 14% for the Antennae). The orbital decay 

6 The eccentricity of 1.10 is not large enough to overcome dy- 
namical friction and leave the galaxies unbound. 
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Fig. 12. — Effect of the eccentricity (top) and the pericenter 
distance (bottom) on the mass fraction in compressive mode. A 
low eccentricity concentrates the compressive events during the 
only passage, while a higher value involves a third passage (at 
~ 975 Myr). A close encounter allows a high mass in compressive 
region. 

brings the two progenitors back at t ~ 700 Myr, and 
the major compressive event occurs. In the close passage 
case (OC), almost the entire mass is available to enter 
into the compressive regime. Therefore, the first peri- 
center peak is much higher (~ 18%). Dynamical friction 
plays a crucial role as it prevents the progenitors to fly 
away from each other. As a consequence, only one long 
compressive episode is visible. 

Here again, we set the end of the integration time to the 
merger stage, to have the same dynamical phase for all 
the orbits. As expected, the short period (ti) does not 
depend on the eccentricity. However, we note a slight 
influence of the pericenter distance (see Table [3]) . The 
distant encounter presents a shorter period (~ 6 Myr in- 
stead of 8 Myr) because the passages are much quicker, 
which is visible with the width of the peaks in Fig. ___J 
For the OC case, the merger stage is too close to the peri- 
center. Therefore, the characteristic timescales derived 
are not physically relevant for a comparison. 

The broad conclusion of this survey is that a change 
of the orbital eccentricity has a mild effect on the mass 
loss and the orbital period. Only the pericenter distance 
influences significantly the importance of the compressive 
modes. To be more specific on this last point, we explore 
in the next subsection the trend of compressive modes 
with impact parameter. 

5.3. Fly-bys versus mergers 

Star clusters also form in galaxies that are undergoing 
milder interaction, as starburst galaxies show different 
levels of intensity. To see how the statistics of compres- 
sive modes is affected by the character of the interaction, 
we varied the impact parameter d of the two progenitors, 
ranging from one to 14 disk scalelengths Rd, keeping all 
other quantities unchanged. In this fashion, the simula- 
tions cover the full range from merging (bound orbits) 
to mildly interacting galaxies (hyperbolic orbits, or fly- 
bys). Nevertheless, our parameter survey remains lim- 
ited in scope, as we do not include high- velocity encoun- 
ters with low impact parameters, which would also pro- 
duce fly-bys. We chose to keep the initial relative veloc- 
ity (~ 200 km.s -1 ) constant, in part because this value 
seems of the right order for field galaxies (as opposed 
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Fig. 13. — Mass fraction in compressive mode at the first pericen- 
ter passage versus the impact parameter d normalized to the disk 
scalelength i?_. Filled circles correspond to compressive modes for 
our runs, while open t riangles are SFR derive d from the hydrody- 
namical simulations of di Matteo ct al. (2007, Fig. 21). 

to ~ 1000 km.s -1 for cluster galaxies); and also because 
higher velocities would lead to shorter timescales, when 
effects of tidal fields may be approp r iately treated in th e 
limit of a shock (see iSpitzerl [19871 iGieies et alJ 120071 ). 
However, let us note for the sake of clarity that a fast 
(shock) encounter will trigger a strong response from the 
gas whenever the impact parameter is so small that the 
disks overlap. 

We graph on Fig. [___ the mass fraction in compres- 
sive mode at the first pericenter passage as a function 
of d/Rd- Each filled circle on the figure represents one 
simulation. The vertical dotted line marks the minimum 
distance for the turnover from merger to fly-by inter- 
actions. Crucially for non-merging galaxies, the mass 
fraction remains at the same level as that of the pro- 
genitors. When d/Rd is decreased, the fraction increases 
slowly but peaks sharply when d/Rd ~ 1. As a result, 
if we were to match compressive modes with star for- 
mation episodes, we would predict enhanced formation 
of stellar associations with decreasing impact parame- 
ter, and no enhancement for fly-bys. It is interesting 
to contrast this wi t h sim ulations including gas. Fig. 21 
of Idi Matteo et all (|2007[ ) shows an enhanced formation 
when compared to fly-bys for all impact parameters in 
the range d/Rd ~ 2 to 6 (Rd is taken equal to the pa- 
rameter a* of their Miyamoto-Nagai models). However, 
these authors also found a trend whereby the star forma- 
tion rate (SFR) increases slightly with impact parame- 
ter, which may be interpreted as a result of more gas 
falling in the deep potential of the final merger phase. 
If we focus on their data for late- type (disk) galaxies, 
this trend is not incompatible with a constant maximum 
SFR, irrespective of whether the orbits are retrograde or 
prograde. This is very akin to our own data for mergers 
in the same parameter range. On Fig. [T__ we display 
as open symbols their data for spiral-spiral mergers (de- 
noted L-L by them). The good overlap of values sug- 
gests compressive modes may be very important during 
a galaxy merger. On the contrary, our data for fly-bys 
indicate no enhancement whereas hydrodynamical runs 
show an enhancement of a factor two or more. This is 
indicative that a fine treatment of the hydrodynamics 
is essential to predict a global SFR, while gravity alone 
would not be enough. Note importantly that we lim- 
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ited our study to scales of 40 pc or larger so that this 
conclusion only applies to relatively massive stellar asso- 
ciations. 

5.4. Mass ratio 

One key parameter in the st udy of interacting pairs 
of galaxies is the ma s s ratio ([Naab k. Burkertl 1 20031 : 
Bournaud et alj l2005t iNaab, Jesseit, fc Burkertl l200a 
Johansson. Naab. &: Burkertl l2009f L Most of the obser- 
vational and theoretical facts set apart major mergers 
(mass ratio greater than 3:1) and minor mergers. To see 
how this could influence the compressive tidal mode, we 
introduce three new models, in addition to the Anten- 
nae (of mass ratio 1:1) reference: another major merger 
(2:1), a minor merger (10:1) and an intermediate mass 
ratio (3:1). 

In this study, the mass ratio refers to the total mass 
(stars + dark matter) and thus, is closer to the c osmolog- 
ical d efinition than to the observational one (see lStewartl 
2009). The galaxies have been scaled to keep the den- 
sity and the virial ratio constant. This choice permits to 
maintain the intrinsic compressive region at ~ 3% of the 
disk mass. 

On the one hand, the small intruders are strongly af- 
fected by their passage in the vicinity of the main disk. 
However the tidal effects they undergo are mainly dis- 
ruptive (extensive) and thus, only their central region 
remains compressive. On the other hand, the main disk 
stays coherent and its mass loss is much less important 
than for its companion. It is important to note that the 
passage of the small galaxy through the large disk in- 
duces radial perturbations which rapidly gives birth to 
spiral arms and then, a bar. These transient features 
generally allows the formation of local cores in the grav- 
itational potential, and thus, the raise of local compres- 
sive fields. 

As a consequence, the mass fraction in compressive 
mode is higher for the main disk, when a small intruder 
hits it, than in isolation. However, the intruder itself 
does not show a strong change in its compressive behav- 
ior. Finally, the total mass fraction in compressive mode 
should be low when small galaxies are involved. 

Indeed, our results (see Fig. [TJ} show that this fraction 
is considerably lower for a small mass ratio (3:1 and 10:1) 
than for the major mergers (1:1 and 2:1). We detected 
features common to all the simulations, such as the pres- 
ence of a peak to denote the creation of the tidal tails, or 
an oscillatory phase for the later steps of the merger. As 
the different masses change the orbital period, the main 
peaks are delayed for the low mass progenitors. 

In the cases where spirals and bars are created, the 
compressive mass fraction presents a plateau during the 
lifetime of these structures. The associated timescales 
are similar to those of the tidal tails and bridges (see 
Fig. [7]). As a result, the exponential timescale t\ (of the 
lus integrated until the same dynamical stage) is slightly 
shorter for the minor mergers than for the major ones 
(see Table El). 

5.5. Further compound models 

The last step of our survey deals with different com- 
pound galaxy models. Here, the main goal is to quantify 
the influence of the halo on the evolution of the merger, 




t [Myr] 

Fig. 14. — Evolution of the mass fraction in compressive mode 
for major (top panel) and minor (bottom panel) mergers. The low 
mass progenitors present much less mass fraction in compressive 
mode than the higher mass galaxies. 

both in terms of the total mass of the halo, i.e. its trun- 
cation radius, and its mass profile. This last point should 
bear on the growth of bridges and the orbits of the nuclei, 
since a more peaked dark matter halo will concentrate 
more mass near the center. It is therefore of interest to 
compare the statistics of compressive modes with mod- 
els where the halo does not have an isothermal core, as 
does the reference Antennae model discussed extensively 
in Section SJ 

At first sight, it should have sufficed to change the halo 
mass profile and keep the stellar content as in the refer- 
ence model to ease comparisons. However, the stability 
of the disk especially requires the coherent adjustment 
of all three components. We made use of equilibrium 
distribution functions from the literature to create new 
composite models, all scaled to have a circular velocity 
of 220 km.s -1 at the solar radius. 

The first model is a Kuijken-Dubin ski model (hereafter 
KD) of a Milky Way-like galaxy flKuiiken fc Dubinskl 
[19951 model MW-A). It uses a lowered Evans distribu- 
tion function for the halo (« 1.3 times more massive 
than the Antennae's), a King model for the bulge and a 
3D version of a Shu disk, truncated at the same radius 
than the reference model, but « 1.5 times more massive. 
As mentioned before, the mass fraction in compressive 
mode in equilibrium should vary from one model to the 
next, as this fraction is linked to the form of the poten- 
tial. The total density in the disk built by the overlap 
of the disk, bulge and halo sets the local Jeans length 
over which a fragmentation mode may grow and form 
cores of compressive tides. Because the Jeans length will 
increase as the density decreases, a more extended mass 
profile, resulting in a shallower potential, is expected to 
yield a lower mass fraction in compressive modes. As the 
KD model is more extended than the reference Antennae 
model, we derived a mass fraction of w 0.8% instead of 
-3%. 

The KD model goes through the by-now usual evolu- 
tionary steps: a first passage, followed by the formation 
of the bridges and tails, before merging completely. Be- 
cause the mass ratio between the two progenitors is the 
same as the reference model, the Keplerian orbits of both 
models match well. Despite the fact that the KD model 
is more extended, differences in the orbits attributable 
to dynamical friction are not significant. Indeed, the 
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pseudo-period of ~ 300 Myr (as denned in Section 15. ip 
is almost the same. However, the low average density 
of the KD model leads to a more effective dispersion of 
the disks material, which results in depopulated nuclei. 
Bearing in mind the considerations of the previous para- 
graph, one may deduce that the compressive mass frac- 
tion will rise, at the first passage, by a factor similar to 
the Antennae model, but by a smaller amount at the 
later stage of merger, because the core region has been 
more effectively depleted (see Table [3]). 

Even if KD and the reference model share a similar 
rotation curve, the energy and angular momentum of in- 
dividual orbits are not the same. We computed shorter 
timescales (~ 2.5 Myr) for the compressive tides, show- 
ing a more efficient orbital mixing in isolation as well 
as during the merger. A detailed explanation for this 
seems out of reach as too many parameters have been 
changed. However, the exponential distribution of the 
compressive mode still has a clear signature, but with 
shorter e-folding times. 

The next model is taken from iMcMillan & Dehnenl 
(|2007l hereafter MD). The halo component is a massive 
NF W profile w i th a s cale length of ~ 12 kpc. The bulge 
is a iHernquistl (|1990l ) profile, while the disk decays ex- 
ponentially with cylindrical radius over a characteristic 
scale w 2.8 kpc. Note that the dark matter halo is much 
more massive and extended than in the previous models 
to the extent that the halos would overlap if set on the 
same orbits used for the Antennae model. Indeed, the 
entire disk sits within on scale length of the NFW profile 
(where p oc Even if this central region is intrin- 

sically tidally extensive, the addition of other potentials 
(disk and bulge) could more easily build compressive ar- 
eas than in the external zones {p oc r~ 3 ). As before, one 
should not expect the intrinsic mass fraction in compres- 
sive mode to be the same as for our reference. In this 
case, the fraction is 0.6%. We noted that not only the 
central part showed compressive tides, but so did a num- 
ber of small areas distributed all over the disk. These 
small zones are clearly compressive because of a tran- 
sient distribution of the matter at small scales, and thus 
do not have long compressive episodes. The same occurs 
in the course of the merger. 

Here again, the time-evolution of the mass fraction 
presents the usual sequence of peaks. We note that the 
high halo mass (24 times that of the disk) reduces con- 
siderably the orbital period. In addition, it increases the 
role of the dynamical friction, as both galaxies stand in- 
side each other's halo. As a consequence, the merger 
stage occurs earlier than for the other models (i.e. ~ 
200 Myr only after the first passage). In the inner part of 
the NFW potential, the overlap of the two disks induces 
many compressive areas. In Fig. |H] (last panel), one ob- 
serves much more compressive particles in the tidal struc- 
tures (bridges and tails). Indeed, between the pericenter 
passages, the compressive regions in the arms do not have 
enough time to dissolve. Consequently, the mass experi- 
encing compressive modes remains fairly high in the MD 
run, while it has time to decrease in the other Antennae 
- and KD simulations. This also imprints the timescales. 
Indeed, the combination of both the transient small com- 
pressive areas and the short orbital period limit the ex- 
ponential timescale of compressive mode to ~ 3.2 Myr. 



The results listed in Table [3] reflects these points. 

In conclusion, even if a modification in the mass dis- 
tribution of the galaxy model does induce changes in 
the statistics of compressive modes, their general behav- 
ior (e.g. location in the tails and bridges, exponential 
timescale, timing) remains overall unchanged, and, espe- 
cially, the statistics of compressive modes at the time of 
the final merger is relatively robust against details of the 
models. 

6. DISCUSSION 

6.1. Statistics of the tidal field 

The gravitationnal tidal field of galaxies is often pre- 
sumed to disrupt smaller bodies orbiting in their vicinity. 
This point of view stems from a wealth of observational 
data of on-going disruption of tidal dwarf galaxies, as 
well as strong theoretical underpinning going all the way 
back to Roche and the notion of an outward stretch oper- 
ating on finite-size bodies. But a close inspection of the 
full tidal field tensor leads one to conclude that in several 
realistic situations, such as in the core of flat-top density 
profiles or the mid-point of spiral disks, the tide will work 
in reverse and act to shelter structures from dissolution. 
Because the character of the tidal field is sensitive to 
the presence of clumps or other sub-structures, it seems 
natural to seek out quantitatively the evolution of the 
tides in the context of a major galaxy mergers, when 
large structures and smallish clumps will form in great 
numbers. A key question that we have addressed in this 
paper is whether compressive tidal modes could widely 
influence the early life of star clusters or TDGs, both in 
terms of their location in the merger, and their duration 
in time. 

The reference Antennae model discussed earlier shows 
that fully compressive tidal modes spread throughout the 
merger, including in key areas such as the nuclei and the 
tidal tails. The statistics of the duration of these modes 
reveal that most of them last long enough to impact on 
the formation of bound substructures on a scale of ~ 
200 pc or smaller. This conclusion was confirmed in a 
parameter survey of galaxy mergers, which showed that 
the shape of the tidal mode distribution function is by far 
unaffected by details of the systems undergoing a merger 
(spin-orbit coupling, mass ratio, impact parameter, and 
so on). 

Broadly speaking, compressive tidal modes always de- 
velop to reach the same level of intensity in all the sim- 
ulations that we have performed. However, some trends 
have been highlighted. First of all, the mass fraction en- 
tering compressive modes scales with the available mass: 
a low-mass progenitor undergoing a collision yields a low 
mass-fraction of compressive tidal modes. This is clear 
from the spin- and mass ratio explorations. Secondly, an 
encounter with a small pericenter distance brings more 
material to the collision in a smaller volume. As a result, 
the number of mass elements experiencing compressive 
tidal modes is significantly higher for close passages. 

6.2. Tidal field and star formation rate 

For every interaction leading to mass exchange and 
possibly a merger, the fraction of bodies experiencing 
a compressive tidal field rises by a factor of at least 4, 
and up to 13 in the extreme case of a face-on collision, 
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when compared with galaxies in isolation. These figures 
match well with the results of Idi Matteo et all (|2008f ) 
who presented a study of the SFR over a large set of 
mergers, using both SPH and sticky particles techniques. 
We also remark that the starburst activity derived from 
their simulations would not last longer than ~ 500 Myr, 
which matches well with our results. This provides in- 
direct evidence of a close link between fluctuations in 
the global galactic gravitational field, and the small-scale 
v olumes where stars a re presumed to form (recall Fig. 2 
of IRenaud et al.ll2008l for another, complementary argu- 
ment^ 

iDi Matteo et all (|2008l ) highlighted that retrograde en- 
counters enhance the SFR of mergers, which is again in 
agreement with the evolution of the tidal field for this 
case (cf. models PR and RR of Section ISTTj) . However, 
they noted a negative correlation between the pericenter 
distance and the starburst, in the sense that close pas- 
sages correspond to low SFR, for mergers only. It is yet 
unclear to what extent this trend can be attributed to 
the response of the gas in the early stages of the merger. 
Because the statistics reported here do not account in 
details for the sequence of events that take place, we can 
not be specific about the likelihood of forming stars at 
a precise time. A clearer picture will emerge once high- 
resolution hydrodynamics is included in the dynamical 
scenarios that we have reported here. 

6.3. Age distribution function 

In all merger simulations that we have carried out, a 
double exponential law describes well the compressive 
mode distribution in age, for both the longest uninter- 
rupted sequence (lus) and the total time (tt). We note 
that the longest periods of the lus statistics correspond to 
orbits near the centre of the galaxies (in isolation as well 
as during the interaction) while the transient features 
that give rise to the shorter periods of the lus distribu- 
tion were found in the merger-induced tidal structures 
(such as bridges and tails). 

In Section [4~51 we showed that the history of the grav- 
itational tidal field along individual orbits is very com- 
plex. Indeed, to compute the net effect of the tides on 
a star cluster over time via a statistical approach, one 
should combine both the lus and tt distributions. This is 
so because the tidal field changes character (compressive 
or extensive) repeatedly along most orbits (see Fig. [5J| . 

To illustrate this further, consider once more the 
case of the Antennae model. Roughly speaking, the 
statistics of compressive modes for that case can be 
split in two phases: a first phase in the time interval 
ATi =]0, 50] Myr, during which the compressive modes 
are mostly uninterrupted; and a second phase AT2 = 
[50, 100] Myr during which the compressive modes are a 
collection of discrete episodes. Putting this together, a 
proxy for the distribution of time spent in compressive 
tidal modes might consist in taking the lus distribution 
over the time interval ATi, followed by the tt statistics 
over the interval AT2. Fig. [15] plots the lus and tt dis- 
tributions (red and blue lines) showing the usual double 
exponential, together with the linear combination of the 
two described above (black curve). 

If we take the bold step of identifying the duration of 
the modes (x-axis on Fig. [TS"]) . with the age of stellar as- 
sociations, then the vertical y-axis may be interpreted as 
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Fig. 15. — The lus and tt distributions (red and blue lines, re- 
spectively) of compressive modes for the Antennae model. The 
distributions consist in double exponential laws fits to the short 
and long timescales (ATi and AT2). The black line represents a 
combination of both cases, using the short timescale from the lus 
and the longer one from tt (see text for details) . The green point s 
with horizontal binning are HST data lifted from Fall ct al. (2005). 

the number of associations per age interval, dN/dr. The 
black curve on Fig. [15] can be fitted with a t~ 108 power 
law (cr = 2 x 10~ 4 ) over the interval t E [10, 100] Myr). 
It is striking that the duration of the compressive modes 
is very s imilar to the r -1 distribution derived from HST 
data by iFall et ail (|2005l ) for the age of young star clus- 
ters. While this does not constitute a full cause-to-effect 
chain of reasoning, the unexpected agreement does offer 
a strong hint that time-dependent tidal fields bear on the 
demographics of young clusters in galaxy mergers. 

6.4. Energy of compressive modes 

To characterize statistically the strength of the tides, 
we plot in Fig [16] the distribution of the maximum eigen- 
value of the tidal tensor for all particles, at selected times. 
As mentioned before, a negative value stands for a fully 
compressive mode. It is clear that, at any stage of the 
merger, the distribution can not be fitted by simple an- 
alytical functions (such as a Gaussian or power laws). 
Its high degree of asymmetry with respect to the origin 
mainly stems from the distribution of eigenvalues in the 
galaxies in isolation (black curve on Fig. ri7j|) . where a 
large fraction of the mass (~ 97%) experiences exten- 
sive but weak tidal modes, while the few percent mass 
elements nearest to the galactic center experience com- 
pressive modes of large, negative eigenvalues. The dis- 
tribution of compressive and extensive modes tends to 
become more even-handed in the course of the merge, in 
terms of the occupation percentage at given |A|. This is 
illustrated at the first pericenter passage, t = Myr, and 
at t = 300 Myr when the merger is in full swing. One can 
deduce from this that the strong destructive effect of the 
extensive tidal field acting at some point in the system 
is balanced out, on average, by a compressive mode of 
comparable strength in some other region, like e.g. the 
tidal tails. Given that individual orbits enter and exit 
compressive modes several times during the merger (cf. 
Fig. [8]) one may anticipate that transitions from com- 
pressive to extensive modes are relatively sharp over a 
timescale of ~ 2 Myr. 

We already noted how the full impact of the compres- 
sive tides could only be quantified once the hydrodynam- 
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Fig. 16. — Normalized distribution (in percent) of the maximum 
eigenvalue (A max ) of the tidal tensor for the Antennae run. Three 
key steps are shown: the configuration in isolation (t = —100 Myr), 
the first pericenter passage (t = Myr) and the merger stage (t = 
300 Myr). 

ics of star formation is included in the modeling. Despite 
this caveat, one can estimate the relative importance of 
compressive tides by comparing the extra binding energy 
they add to a cluster, to the feedback from stellar winds. 
For simplicity, let us assume an isotropic tidal field, i.e. 
described by a tensor set equal to the identity matrix 
times A. For a mass distribution of total mass M t , the 
corresponding binding energy driven by the tide is 
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where a is a dimensionless quantity of the order of unity, 
which encapsulates details of the mass distribution, and 
Rt is the truncation radius. Applying Eq.[3jto the resolu- 
tion of our simulations (M t = 10 M Q and R t = 10 2 pc), 



we get a typical value for A of 



-10 



-30 



(< for a 



compressive tide). This yields an estimate for the tidal 
binding energy of ~ 10 49 ergs. Clearly, this value does 
not compare to the energy release ~ 10 51 ergs of a single 
supernova explosion, and hence SN-driven gas expulsion 
would proceed largely unaffected by the tides. However, 
the binding energy computed from Eq. [3] would corn- 
are favorably to the kinetic energy of an O-s tar's wind 
Cappa fe Benaglial 119981: iMartm et"al] I2008J ). In that 
case the mass loss driven by the winds of early-type stars 
would be slowed down considerably. The enhanced mass 
retention would then help stop or reduce the dissolution 
of young clusters, an aspect that has immediate conse- 
quences to the question related to their mass functions 
in mergers. 

In addition to star clusters, compressive tidal modes 
might also be important for the formation of TDGs. 
Detailed high-resolution ./V-body simulations showed 
that the formation of T DGs is virtually impossible 
in pu rely stellar systems (IWetzstein. Naab. fc Burkertl 
2007). Consequently, a dissipative component (gas) is 
a necessary ingredient to form TDGs. Because the com- 
pressive region will span kiloparsec-size volumes, enough 
stars may form that a TDG will ensue. The binding en- 
ergy of the tidal field estimated from Eq. [3J for a mass 
Mt = 10 9 M Q and a size of 10 kpc (typical values for 
the compressive regions found in the outer parts of our 
merger simulations) is ~ 10 ergs, or the kinetic energy 



release of 10 6 type-II supernovae. Now, if one converts 
10 9 M Q of gas using a standard stellar initial mass func- 
tion and a star formation efficiency (SFE) of 10 — 20%, 
one ends up with on order of a million type-II supernova 
candidate stars. While this calculation is rough, it does 
allow us to conclude that the energy reservoir of com- 
pressive tides must be factored in more detailed analysis 
of TDG formation models, as it is at least comparable to 
SN energy feedback. 

6.5. Perspectives 

In this last, more speculative section we explore briefly 
three key areas for future work that will draw from the 
current study: 

• The first and more obvious is the modeling 
of the time-evolution of cluster mass func- 
tions driven by galactic tides. A number 
of previous studies derived cluster disso- 
lution rates based on steady external tides 
(iBaumgardt fc Makindl2003t IGieles fc Baumgardtl 



2008 : IGieles fc Bastianl 120081: iKupper et alJ 
Vesperini. McMillan, fe Portegies Zwartl 



2008 



2009 



Ernst. Just, fc Spurzeml 120091 ). Other approaches 

time-evolution in an 
the model 



encapsulate the tide's 
approximate "shock" treatme nt as 
cluster orbits the host galaxy (|Spitzerl Il987t see 
IGieles et al.l 120071 for a recent application to spiral 
galaxies). We note that the tidal tensor approach 
taken in this paper encompasses those mentioned 
here as limiting cases. But the important point 
to stress in our view is that the statistics of tidal 
fields change very significantly once we consider 
out-of-equilibrium dynamics. It will be important 
in future to take into account this evolution of the 
tidal field itself when considering the evolution of 
cluster mass distributions, especially in mergers. 

• Secondly, the tidal field bears down on star-forming 
regions at the mass- and length scales of star 
clusters in merging galaxies, and hence on the 
hydrodynamics of fragmentation and star forma- 
tion in such systems. Future modeling of the 
first ~ 10 6 yr of a giant molecular cloud cm- 
bedded in a realistic, fully evolving background 
tidal field should shed new light on questions per- 
taining to th e survival of young (open or globu- 
lar) c lusters (|Bastian fc Goodwin] 120061 : IGoodwir] 
2008). The survival of clusters is set from a com- 
bination of seve ral effects, e.g. the expulsion of 
the residual gas (lGoodwm1ll997l : iGever fc Burkertl 
120011: iChen fc Kol l2008frTand tw o-body relaxation 
(|Portegies Zwart fc Chenl [20081. How the rapid 
evolution of the tides that we have reported 
here changes the picture derived from the sec- 
ular evo lution of clust e rs (the weak field lim it, 
see e.g. iGoodwinl 119971: iBoily fc Kroupal l2003alP 



Bau mgardt fc Kroupal l2007t iParmentier fc Fritza 
2009) is an aspect worthy of further consideration. 



• Finally, several authors have confirmed in re- 
cent years the presence of multiple stellar main 
sequence s (MS) in resolv e d massive c luster s 
(see e.g iPiotto et all 120071: iMilone et ail [2009). 



16 



iD'Antona fc Caloil ()2004[ ) suggested that a stellar 
initial mass function (IMF) with a flat slope at 
the high-mass end would enhance CNO abundances 
yet leave relatively unchanged the abundances of 
iron-peak elements, providing a fit for the multiple 
threads in cluster HR diagrams. It is unclear why 
a bias in the stellar IMF should develop in some 
clusters. One possibility comes in the form of an 
age-spread as broad as ~ 10 Myr, which could be 
interpreted as a second episode of star formation. 
This would open up the possibility that (i) clusters 
are self-enriched (by retaining slow stellar winds); 
or (ii) accrete gas from gas-rich environments met 
on their galactic orbit. In that context it is clear 
that compressive gravitational tides would help in- 
crease the likelihood of a second episode of star 
formation. Gas accretion by the cluster should im- 
pact on both the CNO abundances and those of 
the iron-peak elements, and so offers little hope of 
a viable explanation to this problem. Furthermore, 
the cluster has little time over 10 Myr to move sig- 
nificantly away from the volume where it first con- 
densed before the second burst of star formation 
takes place. Thus, this scenario would preserve the 
homogeneity of the metal abundances the cluster 
already holds. Self-enrichment, on the other hand, 



should be more effective in broadening CNO abun- 
dances since the new abundances would remain cor- 
related with those of the host cluster at birth. Re- 
cent modeling of rotating massive stars has shown 
that slow disk-like ejecta may be more effectively 
retai ned by the cluster than was assumed up t o 
now (jD'Ercole et alj [20081 : iDecressin etHI [20071. 
Clearly the retention of light elements would be en- 
hanced if the cluster was undergoing a compressive 
tidal mode at the time when such winds were oper- 
ative. It is too early to pin down quantitatively the 
actual impact of tidal modes on such cluster evo- 
lutionary processes, which will surely require ex- 
haustive modeling of the hydrodynamics and the 
chemical evolution, not done here. 
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